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,  Abstract 


Tests  of  bomb  stability  require  that  a  large  initial  amplitude  of  yaw 
be  obtained*  The  initial  amplitude  of  yaw  should  he  slightly  greater  than 
45°  but  less  than  90°,  and  accordingly  65°  has  been  considered  desirable* 
The  initial  amplitude  of  yaw  is  dependent  upon  both  the  initial  angle  of 
yaw  and  the  initial  angular  velocity  of  yaw.  Col.  L*  E.  Sinon  has  suggest¬ 
ed  an  apparatus  consisting  of  a  taut  cord  with  one  end  fixed  to  a  pawl  in 
the  bomb  bay  and  with  the  other  end  attached  to  a  ring  engaging  a  cut  nose 
plug.  The  nose  plug  is  cut  in  the  fom  of  an  arc  of  a  circle,  the  arc 
beginning  at  the  junction  of  the  longitudinal  axis  of  the  bomb  and  the  nose 
plug  ring.  Tho  initial  motion  of  the  bomb  in,  and  immediately  below,  the 
bomb  bay  takes  place  in  a  disturbed  airflow.  The  initial  motion  may  be 
regarded,  roughly,'  as  taking  place  in  three  regimes.  In  the  first  regime 
the  air  in  the  bomb  bay  may  be  considered  motionless,  in  the  second  the 
air  in  the  lowest  portion  of  the  bomb  bay  and  just  below  it  as  turbulent 
but  partaking  to  some  extent  of  the  motion  of  the  air  far  below  the  air¬ 
plane,  and  in  the  third  regime  the  air  may  be  considered  as  moving 
horizontally  with  uniform  speed  with  respect  to  the  airplane.  Without  data 
which  are  at  present  unavailable  the  character  of  the  action  of  tho  air  is 
incapable  of  a  precise  description.  As  a  consequence  tho  motion  of  the 
bomb  will  be  considered  as  taking  place  in  air  which  does  *not  move  with 
respect  to  tho  airplane  until  the  retaining  ring  slips  from  the  nose  plug, 
and  thereafter  in  air  having  uniform  horizontal  velocity  with  rospect  to 
the  airplane.  The  effects  of  a  gradual  entry  into  the  free  air  stream 
are  noglected  as  unnecessarily  complicated  for  the  present  problem.  The 
equations  of  motion  for  this  constrained  system  are  worked  out  and  a 
method  found  for  deduction  of  the  point  at  which  the  ring  should  be  cut 
so  that  the  assigned  initial  amplitude  of  yaw  will  be  obtained.  The 
method  is  applied  to  the  Bomb,  Demolition,  1000  lb.,  M44. 


X.  Purpose  of  a  Stability  Teat  for  Bombs 
1*  Aerodynamic  Forces  and  Tor  mes 


The  forces  tending  to  change  the  velocity  of  the 
center  of  mass  of  the  bomb  are  currently  taken  to  be  the  weight 
of  the  bomb,  W,  the  drag,  D,  and  the  lift  or  cross-wind  force, 
L1 •  These  are 

w  =  nv  1 


D  =  Kppdau* 


(1) 


L  =  K^pd2u2  sin  &  J 


where  m  denotes  the  mass  of  the  bomb,  g,  the  acceleration  of 
gravity,  p,  the  density  of  the  air,  d,  the  master  diameter  of 
the  bomb,  u5the  velocity  of  the  center  oX*  riass,  e,  the  angle  of 
yaw,  Kp^the  arag  coefficient  ana  the  lift  coefficient. 

It  will  he  assumed  that  the  boiib  does  not  spin,  tnat  is, 
that  the  homo  does  not  rotate  about  its  longituuina^  axis. 

The  angular  motion  of  the  bono  will  De  sup po sea  to  consist  only 
of  motion  of  the  longitudinal  axis  of  tne  bomb  in  a  vertical 
plane  parallel  to  the  axis  of  the  airplane  at  the  instant  of 
release.  Let  tiie  angle  Detv/een  the  axis  oi  t.,e  bomb  and  the 
horizontal  plane  through  tne  center  of  musa  ue  denoted  by  r\  and 
the  angle  oi*  inclination,  or  angle  between  the  horizontal  plane 
and  the  tangent  to  the  trajectory  be  denoted  by  9.  It  appears 

from  Diagram  1  that 


DIAGRAM  1 


T}  =  £  +  0  •  (2) 

Let  T  denote  tne  torque 
acting  to  turn  the  bomb 
about  its  center  of  mass 
and  I.  the  mo-.ent  of 
inertia  of  tne  bomb  about 
a  transverse  axis  through 
the  center  of  mass,  by 
the  buierian  physical 
equation2 : 


l 

The  formulae  for  forces  and  torques  given  in  this  report  are,  with  two 
exceptions,  those  employed  in  Ballistic  Research  Laboratory 

Report  No.  Ill  "Aerodynamic  Nomenclature  and  Formulas,  Conversion  Factors  and 
Drag  Functions".  The  angle  of  plane  yaw  is  positive  when  the  nose  of  the 
bomb  is  belcv,  the  tangent  to  the  trajectory,  and  the  expression  for  lift  has 
the  sign  of  the  angle  of  yaw.  The  damping  torque  denoted  by  Hitchcock  as  Hco 
is  here  denoted  by  H  in  conformity  with  the  well-known ‘work  of  R.  H.  Fowler. 
It  is  now  considered  that  the  force  system  here  given  is  somewhat  incomplete 
and  a  thorough  study  of  this  system  is  now,  being  made  by  J,  Kctihane. 
a 

J.  H.  Jeans  in  "Theoretical  Mechanics"  derives  Euler* s  physical  equation  on 
pages  306-303  and  pages  346-348, 


--JK- 


T  6ft  _  T  # 

IB  dt  ~  T 


The  angular  velocity  of  the  axis  of  the  bomb  with  respect  to 
directions  fixed  7/ith  respect  to  the  air  is  denoted  by  -  Q  . 
Evidently 

C-T)  =  e  +  0*  (4) 

The  external  torque,  t  ,  is  the  sum  of  two  torques,  the 
restoring  torque,  denoted  by  M  and  the  damping  torque  denoted 
by  H.  These  are 

M  =  -  Kjjjpd^u*  sin  e  1 


H  =  - 


Kj^pd^uQ 


Substitution  of  the  value  of  t  in  the  Eulerian  Physical  Equation 
(3)  results  in 

IB  {f-  =  -  KMpd3u2  sin  -  KHpd4uC  •  (6) 

Substitution  of  the  value  of  £  from  equation  (4)  in  equation  (6), 
transposition  of  terms  ana  division  of  both  sides  of  the  result¬ 
ing  equation  by  yields 

•.  KMpd*u2  Kupd^u(fc+0) 

e  =  -  - sin  e  -  — — z -  -9  •  (7) 

XB 


Let  x*  denote  file  air  range  coordinate  of  uhe  center  of 
mass  counted  forward  from  the  release  point  and  y1  the 
corresponding  depth.  Unaer  standard  bombing  table  conditions 
the  equations  of  motion  of  the  center  of  mass  are  given  by1 s 

**,  _  D  cos  Q  L  sin  9  ^ 

“  "  ra  )  (8) 


t  -  _ 


D  sin  9  .  L  cos  9  . 
m  m  6 


The  angle  of.  inclination  of  the  tangent  to  the  trajectory  has 
been  denoted  by  0.  Then: 

9  =  arc  tan  (£)  •  (9) 

The  equation  (7)^  which  governs  the  motion  of  the  axis  of 
the  bomb  in  the  vertical  plane,  may  be  simplified  considerably 
under  the  conditions  prevailing  during  the  short  time  interval 
which  is  of  present  concern.  The  time  interval  during  which 
the  motion  of  the  axis  of  the  bomb  is  of  interest  covers  the 
time  between  release  of  the  bomb  nose  plug  from  the  retaining 
ring  and  the  time  when  the  bomb  attains  its  first  maximum  yaw, 
or  initial  amplitude.  This  interval  is  iess  than  a  quarter 

1  These  equations  are  given,  for  example,  by  F.  V.  Reno  in  Ballistic 
Research  Laboratory  Report  No.  145  "On  the  Theory  of  Motion  of  the  Bomb". 
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period  of  the  yaw,  and  within  this  interval  the  center  of  mass 
makes  a  vertical  drop  of  less  than  six  feet.  The  shortness  of 
this  interval  makes  it  possible  to  neglect  the  variation  of  the 
density  of  the  air,  p.  Accordingly  p  will  be  replaced  by  pQ 
in  equation  (7). 

It  is  found  readily  that  '  N 

KDp0d*u* 

u  =  g  sin  0  -  -  / 


A  =  g  cos  e  +  Vod’u.sl"  6 
w  u  m 

..  r  .  ?rL  ?kl 

0  =[-  KLg  sin  6  +  uu  ~  +  ue  —  + 
epod»uKLc0S  e  '  g,  sln  q  cos  e 
in  "  u2 


1  P«d* 

»]  ^ 


cos  0 
—  u 


sin  e 


The  result  of  substitution  of  (10)  into  (7)  is: 


V  ,  ,  1  ?kt-  , 

mIB  mu 


sin  e  ?kl 
mu  au 


Lp/u?K.  Kt  Kn  . 

-  -V-  ^  “  -5*®  >o<4d^ 

[K,  KHd2_| 

a  cos  e  +  -j-j  P0 d2ue 

.[5^1.^'.2^]6Cos6 


g  cos  0  • 


Mean  aerodynamic  coefficients,  considered  independent  of 

the  an&le  of  yaw  and  the  air  speed,  will  be  employed  in  the 

present  report#  Consequently  aKr  dXT  will  be  considered 

— rif  and 
06  du 

zero,  sin  e  will  be  taken  equal  to  e  and  cos  e  =  1.  The  result 
is*  i-  i  — i 


.  .  fefi  .  KhW4  .  ilWll  .  4, 

L1*  uIB  ~  J  0 

■  fr* +  ^17]  w2** 

}~ KHpoq4  2g  sin  0  KPpou*l 

"L  XB  “  u*  n  J6 


+  goose. 


I 


nwudkt.  «  JL«C 
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For  the  short  period  in  time  with  which  the  present  problem 
is  concerned,  "g8  si^gf  t;0—  is  negligible,  cos  9  may  be 
replaced  by  unity  and  the  airspeed,  u,  replaced  by  uQ  without 
appreciable  error1.  The  result  is;  > 

.  _  rv .  %wl  ......  i 


»•  _ 
e  ®  - 


+  IB. 


p.d2u*e 

rO  0 


KuPJi4  KnPrtd 


.  *HKo  .  Dho 
p/u£  -  -  v  '  +  —  g 

ro  o  ig  m 


The  moment  of  inertia,  Ig  will  be  replaced  by  its  equivalent  in 
terms  of  mass  and  radius  of  gyration,  mk£.  Brief  numerical  com¬ 
putation  indicates  that  KHKLpod  and  KP  vo  are  negligible  in 


comparison  to  and  these  terms  may  be  dropped.  The  result  iss 
B  r-  — i 


P«dX  K, 


e  =  - 


M2  P„d2», 


(14) 


Kj,d 


«P0a‘ 


A  determination  of  the  degree  of  stability  of  flight  can  be 
made  for  any  given  set  of  conditions  if  the  aerodynamic  coefficients 
are  known  as  functions  of  the  yaw  and  Mach  number.  The  purpose  of 
a  stability  test  is  not  to  test  directly  the  stability  of  the  bomb, 
but  to  measure  the  aerodynamic  coefficients.  In  performing  a 
stability  test,  it  is  desirable  to  have  an  initial  amplitude  of  yaw, 
designated  hereafter  as  A,  of  at  least  30  times  the  probable  error 
of  a  measurement  of  the  yaw.  On  the  other  hand,  the  yaw  should  not 
exceed  90°  as  an  overturning  bomb  requires  an  inconvenient  procedure 
in  reduction  of  angular  data.  The  most  suitable  initial  amplitude 
of  yaw  is  between  45°  and  90°,  preferably  about  65°.  It  is  proposed, 
in  order  to  produce  an  initial  amplitude  of  yaw  of  any  desired 
magnitude,  that  the  bomb  be  suspended  from  the  bomb  bay  by  a  cord 


1  No  considerable  approximation  in  the  equation  for  plane  yawing  has  been  made 
previous  to  ll3J,  at  least  for  small  angles  of  yaw.  The  solution  of  (13)  is  by 
no  means  sufficiently  accurate  for  problems  other  than  that  examined  in  this 
report.  A  solution  of  a  more  accurate  character  has  been  made  by  E.  J.  Mc&hane. 
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ring  attached  to  the  nose  oi'  the  bomb«  The  nose  plug  ring  will 
be  cut  at  an  angle  o>  from  tin  axis  of  the  bomb,  as  shown  In 
Diagram  (2)(b).  The  angle  cowill  be  chosen  such  t'^t  the  initial 
angular  velocity  of  yav:  and  iiitial  angle  of  yaw;  will  combine  to 
produce  the  desired  initial  amplitude  of  yaw*  \ 


Eandaixi  Bobbin rt  Table 


Demolition  bombs  are  suspended  horizontally  in  the  bomb  bay. 
As  the  axis  of  the  bomb  is  in  the  same  direction  as  the  line  of 
flight,  the  initial  angle  ox  yaw  of  the  bomb  is  zero.  At  the 
instant  of  release  Q  is  also  zero,  and  the  time  rate  of  change 
of  e  is  equal  to  tnat  of  9,  but  opposite  in  sense.  The  normal 
initial  conditions  for  the  horizontal  suspension  of  demolition 
bombs  ares 

•  ■  °  i 

‘.■-v  j  (i5> 

rfith  regard  to  stability  tests,  the  conditions  are  more  general. 
The  origin  of  time  will  be  taken  as  the  instant  the  bomb  is 
released  from  the  retaining  arc.  The  angle  of  yaw,  e,  the  angle 
of  inclination  of  the  tangent  to  the  trajectory,  0,  ana  the 
angle  of  inclination  of  the  axis  of  the  bomb  are  related  by 

7)  =  e  +  0  • 

Between  the  instant  of  release  from  the  shackle  and  the  instant 
of  release  from  the  retaining  arc,  the  angle  7)  is  the  negative 
of  the  angle  cp.  Then 


S  =  T)-9  =  -<p-9 
4  =  7)-G  =  9-9 


Under  the  assumption  that  the  air 
forces  are  negligible  during  this 
regime,  it  follows  that 

y  =i(cos\|/-  cosv_T)  +  h  sin  <p  i  (17) 
x  =  uQt  -  l(sin^-  slny)  +  h(l  j  co s<? 

Evaluation  of  y  .  k  ,  y  .  x  .  0 
.  •'o'  o'  o'  o  o 

and  60can  be  made  easily  from  (17). 


DIAGRAM  3 
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The  results  are: 


y0  =  l*pQ  cos  <p0  -  h  0  sin 

*o  ='  uo  +  h^o  sin  1  i  003  \ 


••  •• 


y0  =  cos  sin  v0  “  l\  sin  cos  *£> 

'%>  =  h<i>o  sin  %  +  h,?0  cos  %*  cos  %  -  i  sin  ^ 

a  _ .  ..  h%  003  Vo  - 1  i  sin  k  (  (ii 


9  =  arc  tan 


u  +  h6  sin  G  +  IT  cos 


uoh(0ocos  (?Q  -  (p2  sin  <?0)-  u0Z(i*o  sin  y*2  cos  »o)-h20o3 

•  lc0s  (°Q- U  [?''oV'^U-hlsiU  ^c,-Q  [6ag*^»o1 

O  ^o^o+^'V'o  sin  (%"  ^+uo+2Vh0o  sln  Vj*o  cos  V 

Tlie  angle  Q  is,  for  present  purposes,  small.  Brief  numerical 
calculation,  employing  values  given  later  in  this  report,  shows  that 
all  _  in  Qq  or  9q  which  do  not  contain  uQ  are  negligible  in 

comparison  with  terms  that  contain  uQ* Then 

hi  cos  o  -  2i  sin  ^  / 

eo  ‘•’o  u.  ( 


6  =  -  <?  - 
O  Yo 


h(<PQ  cos  <?Q  ~  92  sin  90)-UV0  sin  y  %  oos 

u°  -S 


)T,oxi.T,ate  Deters 


'ej’ticai  Jjiai 


of  the 


J  ■»-*.!  ■  ,  U  * 


-  ■  K  *  f:  1  ^ 


afl.  „  ..  „  ,  %d‘  6|)oa‘ 
H  »  p  "  [S  k£  J  Hi 


Substitution  of  (20)  in  (14)  yields 

e  +  2su  6  +  n*u*(e  +  zrfrir)  =  0 


(21) 


Write  z  for  e  +  ■  and  (21)  becomes 


z  ♦  2suqz  +  n2u2z  =  0  .  (22) 

This  is  the  well-known  equation  for  free  vibrations  of  u 
damped  oscillator.  The  solution  Is:1 


z  —  Ae 


-su  t  r  _ _ 

0  cos  v*i2  -  s2 


uot  -  ®]  ‘ 


where 


,  /z2  +  2u  sz  z  +  z2u2n2 

A_  1  /  O _ OOP  O  O 

"  ..  1/  _  ?  .3  * 


n2  -  s2 


An  easy  numerical  calculation  shows  that  is  of  the  order  of 

n  u* 

0  •  O  e 

0.1  anu  therefore  negligible,  anu  zQ  is  of  course  equal  to  eQ. 
The  solution  for  (22)  is: 


-su  t 

e  =  Ae  *  °  cos 


[y'n2  -  s2  uQt  -  o} 


where 


i  Ai  +  2u  se,£  +  e2u2n2 
a  -  I  o  -  Q  Q  o  o  o 

A  ~  u  V  n2  -  s2 


a  =  arc  sin 


(  £o  +  uoS£o 

\jrrrzrzrr' 


£2  +  2u  +  e2u2n2 

o  o  o  o  oo 


file  natural  frequency  is; 
f  =  u0*/ n2  -  s2  * 


The  period  is: 

P  =  Q  =  20—  —  •  (29) 

uQy  a2  -  s2 

It  is  obvious  from  Diagram  4  that: 

A,  *  Ae  *  SV’'  . 

and  .  _  .  -  su  t,  f  (3°) 

ii2  —  au  O  *  1 

1  A  discussion  of  this  well-known  equation  and  its  solution  nay  be  found  in 
"Introduction  to  Theoretical  Physics"  by  Leigh  Page,  D.  Van  Nostrand  Co.  Inc., 
1928  on  page  63  ♦t  seq. 


(31) 


Solving  equations  (30)  for  s  yields: 


1 

**0 


l°g  6  ^  * 


Now: 


and  therefore 

»■  -  s>  .  (j=-)* 

o 

n*  =  (fjH*  ♦  s*  * 

o 


(32) 

(33) 


The  true  curve  of  yaw  against  time 
will  be  approximately  represented 
by  the  solution  (25)  for  a  very 
short  time,  about  one  second*  As 
the  vibration  is  camped,  the  second 
amplitude  of  yaw,  A2,  is  less  than 
the  first  ,A, • 


DIAGRAM  4 

The  stability  test  yields  observations  of  the  period  and 
the  amplitude.  Measures  of  the  period  anc  the  amplitude  can  be 
used  to  derive  the  auxiliaries  n2  and  s.  The  aerodynamic 
coefficients  ana  +  K^k^2  can  then  be  found  from  the  formulae: 


d* 


A  graphical  procedure  is  convenient  for  the  derivation  of 
A,,  A2  and  P.  The  quantities  directly  measurable  from  the 
photographs  are  the  depth, y,  the  trail,  X  ,  and  the  angles 9  as 
functions  of  count  number,  c.  It  v/ill  be  supposed  that  the  time 
bears  a  linear  relation  to  the  count  number.  The  constants  :>f 
the  linear  relation  may  be  found  by  a  determination  of  the 
origin  of  tine  and  a  determination  of  the  rate  of  the  camera, 
the  latter  being  found  by  photographing  the  face  of  a  stop-watch. 
Then  a  plot  gives: 


y  =  y(c) 
x  =x  (c) 
n  =»)  (c) 


(35) 
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A  plot  of  0  Is  required.  Then 


0  =  arc  tan  ^  =  arc  tan  ^ 


=  arc  tan 


UQAt-AX 


(36) 


The  increments  Ay,  At  and  AX  having  been  read  off  the  plot,  0 
can  be  calculated  as  a  function  of  count  number.  Reference  to 


A,  =  >1,  -  0( 
A2  =  ^2  “  ^2 


] 


(37) 


DIAGRAM  5 


The  density  of  the  air  at  the  altitude 
of  release  may  be  deduced  from: 

p  .  ^  ■ 

where  p  is  given  in  pounds  per  cubic 
foot  and  Y  is  the  altitude  of  release 
in  feet. 


In  the  present  study  the  amplitude  of  yaw  for  the  stability 
test  will  be  assigned.  This  is  done  by  substitution  in  equation 
(26):  _ _ _ 

A  =v  e  2  +  2u  se  c,  +  e2u2n2 
»  o _ o  o  o  o  o 


UqV^ n2  -  s2 

II.  The  Bomb  Constrained  by  a  Cord  From  a  Point  Fixed  in  The 
Airplane  to  a  Cut  hose  Plu£ 

1.  Geometry  of  Constraint  Imposed  on  the  Bomb 


Diagram  6(a)  is  a  schematic  representation  of  the  constraint 
to  be  imposed  on  the  bomb*  At  B,  there  is  a  cut  ring  attached  to  the 
nose  of  the  bomb.  On  this  ring,  will  be  placed  another  ring  with 
a  cord  attached  thereto.  This  cord,  OB,  is  attached  to  a  fixed 
point  in  the  bomb  bay  of  the  airplane  as  shown  in  Diagram  6(b)* 

For  purposes  of  the  subsequent  development,  it  is  assumed 
that  the  cord  is  inextensible  and  acts  as  a  rigid  body  similar 
to  a  very  thin  bar. 

There  is  a  rotation  of  OB  about  0  and  also  a  rotation  of  the 
bomb  about  B.  This  latter  rotation  might  also  be  considered  as 
the  rotation  oi  the  bomb  about  its  center  of  mass.  The  problem 
to  be  solvec  is  thus  one  of  a  system  of  constrained  motion  having 
two  degrees  of  freedom,  both  of  which  are  rotational. 

Let  the  angle  Detween  the  cord  and  the  vertical  be  denoted 
by  y  and  tne  distance  of  tiie  center  jo f  mass  from  the  nose  of  the 
bomb  by  h  .  net  «*>  be  the  angle  between  the  axis  of  the  bomb  and 
the  horizontal. 


Che  Equations  of  Motion  and  Their  Solutior 


fo-rfl  f  t 
u i 


/(/-  *+*9) 


A  most  direct  and  comparative! 
simple  method  of  developing .the 
differential  equations  01  motion 
for  this  system  can  De  found  by 
use  .of  Lagranges*  equations.1 
The  kinetic  energy  of  the  system 
is  partly  due  to  translational 
motion  and  partly  due  to 
rotational  motion.  The  total 

r  kinetic  energy  of  the  bomb  is 
given  by:2 

T  =|m[p*  ♦  q*]  +|lB«p.  (38) 

(c)  L 


The  potential  energy,  V,  consists  of  two  terns,  the  sum  of 
these  terns  being  due  to  tne  height  of  the  center  of  muss  of  the 
bomb  above  its  lowest  possible  position.  The  two  terms  are  due 
to  the  two  distances  l  (1  -  cos y  )  and  h(l  -  sin  p)  which  nake  up 
the  co  ibined  neigivt  of  the  center  of  mass  above  its  lowest  possible 
position.  t..s  ci stances  are  shown  oy  Diagram.  7(b)  and.  Diagram 

7(c).  evidently: 

r  — r , 

V  =  mg  [2(1  -  cos  y)  +  h(l  -  sin  p)J  *  (39) 

1  An  extended  discussion  of  Lagrange's  Equations  is  found  in  E.  J.  Routh, 
"Dynamics  of  a  System  of  Rigid  Bodies",  Part  I,  Chapter  VIII. 

The  coordinate,  p,  of  equations  (38)  et  seq.,  is  unrelated  to  the  temporary 
parameter,  p,  introduced  in  Section  I,  3# 
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By  definition: 

IB  =  “kl  (40) 

vhere  kg  is  the  radius  of  gyration  of  the  bomb  about  a  transverse 

axis  through  the  center  of  mass*  The  Lagrangian  function,  L,  is 
the  difference  between  the  kinetic  and  potential  energies. 

i .  *  - » = a  ■  [p> .  ;*] .  §  iB 

-  mg  [id  -  cos  y  )  +  h(l  -  sin  <?)j  *  (41) 

It  is  obvious  fron  Figure  7(a)  that: 
p  =  -  l  sin  ij,  +  h  cos  cp 

q  =  +  l cos  y  +  h  sin  (?  .  [  (42) 


Therefore: 


p  -  -  cos  v  -  hep  sin  cp  j 

«>  .  •  /  ,  x 

q  -  -  m  sin  y  +  ht?  cos  cp  ,  I  (43) 

The  result  of  squaring  and  adding  equations  (43)  is: 

*  •  >  •  «  •  • 

p2  +  q2  =£y2  +  h2cp2  -  2hZcp^  (siny  cosq>-  cos^  sin  cp) 

•  •  • •  (44) 

=  2%  2  +  h2u>2  “  2hic?4"  sin  -<p)  • 

Substitution  of  this  value  in  equation  (41)  yields: 

L  =  ^  raj?2'*'2  +  h2(p2  -  2h£q>¥  sin  (i'  -cp)] 

+  2  IB<p2  ~  rag  [j^1  ~  cosy)  +  h^1  *  sin  (45! 

=  2  m  02  +  h2)^a  “  2hl(p  ^  sin  (v -<p)j 

+  mgh  sin  <p  +  mg 2  cos  y-  mgh  -  mg 2  . 

• ,  /i* . .  - 

Lagrange^  _  equations  for  finite  forces  are:  1 

-d  (2li)  =  3±  •  1 

W  f  (46) 

and  if  ^  =lr  *  J 

The  required  derivatives  for  Lagrange’s  equations  are  found  readily. 
After  substitution  of  these  derivatives  into  equations  (46)  and 
slight  reduction  it  is  found  that: 

1  Footnote,  P.  11,  supra. 
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n  [k2  +  h2]  <p  =  mgh  cos  <p  +  rah2*  sin  (y-<p)  ) 

L  0  I  (47) 

+  mh cya  cos  (y-<p)  f 

m «  -  mg  2  sin*  ♦  mh2  <p  sin  (*m?)  I 

-  mh2<p*  cos  (v— q>)  •  " 

As  a  check  upon  these  results  it  may  be  noted  that  if  l  =  0, 
the  first  of  equations  (47)  becomes: 

m  |k*  +  h2J  <p  *  mgh  cos  <p  , 

which  is  the  equation  for  a  compound  pendulum.  Also  it  may  be 
noted  that  if  h  =  0,  the  second  of  equations  (47)  becomes: 

mi2Y  =  -  mg  l  sin  V, 

which  is  the  equation  for  a  simple  pendulum. 

A  slight  modification  of  equation;*  (47)  will  be  advantageous. 
Evidently:  _ 

M  •*  y,  /  •  f 

ip  =  ■*“ -  cos  ip  +  fr® —  Y>  sin  (<mp)  +  “ -  ¥*.  cos  ( 


k2+h2 


k2+h2 


k2+h2 


^  sin  V+  y  <p  sin  (V-<p)  -  y  ip2  cos  (y-<p) 
The  angles  <|>  and  y  are  small  in  the  sense  that: 


;os  (h^<p) 


(48) 


•  „  2 

sin  (p  =  <p  cos  <p  =  1  -  £7- 

•  *2  (49) 

sin  v  =  v  cos  y  =  1  -  jjj-  • 

Equations (48)  become: 

<p  -  f& ft. -  -  -  cos  <p  +  ^ ^  ■  (if,  COS  P  —  CD  COS  Uj)  | 

k2+h2  k  2+h2  I 

+  “ —  l|> 2  (cos  If/  COS  <p  +  lf/Cp)  l  (50) 

k2+h2  7 

y  s  -  y  +  J  <P  (y  COS  P  —  p  COS  y)  -  y  92(coS  ty  COS  <p  +  i|/(p) 

*  *  *  J 

The  equations  (50)  are  the  differential  equation®  of JimotionLOf 
thevsystem  unaer  investigation. 


(49) 
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The  initial  conditions  aret 
■  • 

9(0)  =  <?0  =  0 
<p(0)  =  90  -  o 


f  (0)  «  ^ 

*  (°)  =  »0  -  0  . 


(51) 


solution  for  these  equations  will  nov/  be 
developed.  Assume  the  following  power  series  fora  for  q>  and y  s 


»  5  ?0  +  V  +  a»t2  +  a5t3  +  •• 

*  s  +  V  +  b2t*  +  b,ts  +  ..  . 


(52) 


or  +hon?«??^i0r*‘'  11137  be  simplified,  slightly  by  substitution 

are:h  iniUal  conQitlons  as  «iven  by  equations  (51).  The  results 


<P  * 


V  =  YL  +  b2t4  +  b.t3  +  b.t^  +  .. 

Q  +  ^  1  •  •  f 


(53) 


For  the  small  interval  involved  here,  approximately  one 
quarter  of  a  secona,  both  series  are  convergent?  Sver°  as 

fourth^rder  °*2*  secona>  terms  of  higher  than  the 

0r??r  must  be  taKen  mto  account,  in  oraer  to  obtain  a 
sufficiently  accurate  result  for  present  purposes*  On  the  other 

hana,  fourth  order  terras  may  be  neglected  for  very  small  values 
ox  t« 


Differentiation "of  equations  (53)  leads  to: 
9  =  2a2t  +  3a  3t*  +  4a^tJ  +  ... 

*  =  2b2t  +  3b3t*  +  4b4t3  +  ... 


(54) 


and 


•• 

9  =  2a2  +  6a3t  +  12a, ta  + 
••  ^ 

=  2ba  +  6b.t  +  12b. t*  + 

4 


•  1  • 


•  •  •  # 


(55) 


(53h  (5^)  and  (55)  in  the. differential • 
equations  (5u)  ana  neglecting  ail  terns  ol'; 
or  higher  oraers  in  the  time  results  in:  . 


/' 


-14- 


I 


•  •  • 


2a 2  +  6a jt  +  12a, t*  + 

4 

=  +  fa*  V6b>  »ot+12b4  Vi+2b*  [ba  -aJljfil 

+  ...  ■*- 

2ba  +  6b 3t  +  I2b4t*  +  ...  (56) 

=  ¥  fr0+b.f)*J  |&.^,*0t+12a  fot*+2aa(ba-a1flJfe]J 

-f  [^(1.^)]  +  ...  1  2ji 

of  a  function  LTJoKr^eSlHs  Si?*  Sta£es  tliat  the  development 
of  tae  sane  powers  of  t  on  both  sid^o^ecu^??^0??^1*  coefficients 
the  sane.  Equating  coefficients  o-^eqtins  ^  must  be 
slues  of  (56)  leads  to  g™  £d*C 

2a*  =  kW  +  jib,  {| 

6a*  =  —  FT  -l  ..  ""1 


”“J  "  P+P  £mo] 

^4  ■  Mh^v.*2t’[^a-4jJ*  bkrU»ia-*‘)] 

2b*  =  -  f  *0  +  f(2aa»0) 

6b3  =  (6a3»o) 

12b4  =  f  b,  .  ^[l2a,(()  *  2llj  (bj  .  a  f,  .  0,1 

-r[«i  u-Y)]-  1  SJ 

The  latter  is  the  case  here.  ca  3  9  * 


Now 


a*  =  2?k*+h*)  +  P7P  (h*  *0] 
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■  ,  ./'•"V 


b* =  -  h  % +  r  * 


Simultaneous  solution  or  these  two  equations  results  In: 

gh  l-  ^  ^ 

a*  s  / 

"6  1)kS  4  f 

ba  *  22[k2+h2-h2  t£]  *  \ 


Similarly: 


a4  =  LV'+h*)  [12Vo  +  2b2  [b*  -  32  [ 

12(k2+h2)  [fb2  ^  "  “til 


1  -  i2 
1  2 


|153> 


‘4  =  fh  "•  *  fez  [12Vo  *  2<“  [b=  -  “‘l1  -  If]] 

-  fe?  [4»»(i  -  ^j)]  •  J 

The  solution  1  of  these  equations  is: 

a4  =  2gi>hk^  V  (4- v]-3^g2h^»0  {2-5  ^g2haicg  s»Q  [2-3S'0a-^  Sn^j 

4ol  Qc2+h2-h2  yQ*J  3  ^ 

2g2k4  ^(k2+h2)+s2h2k  V2{2k2  ^(3-1-2.)}  -3h2^  >2-5  ^nV4-^! 

b4=  -■  4  ^2^+h^.,%]3 

.  '  Wk2  16  -  17 »Q2  -  15* n4  -  **QUJ  ...■“  - . .  . . 

48 2i  ]F  +  h*  -  h2^’^  • 

®  1 

The  leading  terras  of  the  initial  motion  solution  are  obtained 
by  substitution  of  these  values  for  a2,  a  ,  b2  and  b,  in  the 
following  equations:  4  4 

<?  =  a2t2  +  a^t^  -  . .  •  ^ 

V  =  v0  +  b*t2  +  b^t4  +  ...  I  ^ 

In  uddition  to  the  foregoing  solution  by  power  series 
development  in  the  indepcaueit  variable,  tnere  are  other  possible 
methods  of  solution.  But,  iiasnuch  as  only  tiro  initial  motion  is 
of  interest  for  t.ie  present  p-oblem,  more  powerful  and  general 
methods  of  wOivlnj  differential  e-.uati'jnj  v/ill  not  be  enwijyea, 

’  Tha  viMon  rf  •>,  avl  I*  arc  fir*  r  In  (ss).  Then  tho  values  of 

a^  and  b^  are  obtained  by  slmltarnus  solution. 
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The  equations  of  motion  may  also  be  solved  by  numerical 
integration  in  any  particular  case.  The  method  of  numerical 
integration  to  be  employed  in  an  example  later  in  this  report 
is  referred  to  as  Moulton»s  Method.1  The  result  of  solution 
by  numerical  integration  and  the  initial  lines  of  the  numerical 
work  will  be  given  later.  The  example  given  is  that  of  the 
BOMB,  Demolition,  1000  lb.,  M44*  hut  the  methoa  is  general  and 
can  oe  applied  to  bombs  other  than  the  M44* 

III.  Application  to  the  BOMB.  Demolition.  1000  lb..  M4A 

The  proposed  physical  plan  with  respect  to  the  BdMB, 
Demolition,  1000  lb.,  M44  is  shown  in  Diagram  8. 


The  equations  of  motion  for  this  bomb  are: 

W  M 

<p  -  10.9375cos<p  +  1.0040 y(ycos<p-<pcos y) 

.  M 

Y  =  -10.8942  Y  +  0.7499<p(ycos<?-<pcos  - 

The  computing  schedule  for  numerical  integration  and  several 
lines  of  computation  are  shown  in  Table  I.  These  can  be  used 
as  models  in  any  subsequent  investigations. 

The  numerical  integration  was  continued  to  a  time  equal  to 
about  one  third  of  a  second  to  make  certain  tnat  it  covered  a 
range  resulting  in  values  of  A  exceeding  65°. 

1  The  proof  of  this  method  was  first  worked  out  by  P.  B.  Moulton  in  1914* 
The  method  and  proof  were  published  in  1913  in  a  blue-print  pamphlet  issued 
Vy  the  Ballistic  Branch  of  the  Ordnance  Dept.,  U.S.A.  and  were  reproduced 
in  1926  in  Moulton's  "New  Method  in  Exterior  Ballistics",  Chapter  V. 


+  1.0040^  (cos  Hcoscp+'J* 
0.7499<p2  (cosi|/Cosq)+\fp^ 


|*y?('l'coB'|Mpco8  4  -jp**  vcos  yco8qH-  wp) 


t 

(sec.) 


T 

(rad.  /  sec?) 


4- 

(rad.  /  sec?) 


(rad?  /  sec?)  (rad.  /  sec?) 


0.00 

-2.483 

1.844 

-0.000 

-  0.64 

0.01 

-2.483 

1.840 

-0.009 

-  0.65 

-1 

0.02 

-2.482 

1.825 

-0.034 

-:0.69 

-4 

-3 

By  substitution  of  the  numerical  coefficients  in  equations 
(53)  ^  (54)  and  (55) ,  the  followin z  results  are  obtained  for  the 
initial  motion  solution: 


9  *  5.3957ta  -  2.1700^.+  ... 

9  =  0.2279  -  0.3192ta  -  10.93C0t4  + 


,  (62) 
(rad.) 


9  *  2  x  5.3957t  -  4  x  2.1700t>  +  . 
9  =  -2  x  0.3192t  -  l  x  10.9300t3  + 


.  (rad. /sec.  f) 
. . . (rad. /sec.) 


(63) 


9  =  2  x  5.3957  -  12  x  2.1700ta  +  ...  (rad./sec.a) 
9=  -2  x  0.3192  -  12  x  10.9300ta  +  ...(rad./sec.a)/ 


(64) 


Table  II  gives  a  comparison  of  the  results  obtained  by 
numerical  integration  anu  by  substitution  in  the  analytical 
solution  at  one  tentn  and  at  two  tenths  01  a  second* 


tabu;  i: 

r 

X 

Numerical  Integration: 

t 

(sec  ) 

<P 

(rad) 

$ 

(rad/ sec) 

<!>' 

(rad/seca) 

¥ 

(rau) 

(rad/ sec) 

•  • 

(raa/sec i) 

o.io 

0.-0537 

1.071 

10.60 

0.2237 

-0.107 

-1.92 

0.20 

0.2138 

2.137 

11.00 

0.1981 

-0.468 

-5.71 

Analytical  Solution: 

t 

(sec) 

<P 

(rad) 

$ 

(rac/sec) 

(rad/sec?) 

V 

(rad) 

• 

(rad/sec) 

#  # 

Y 

(rad/seci) 

0.10 

0.0537 

1.070+ 

10.53 

0.2236 

-0.108 

-1.95 

0.20 

0.2124 

2.089 

9.75 

0.1976 

-0.477 

-5.88 

It  will  be  noted  that  at  a  time  equal  to  one  tenth  of  a 
second,  the  two  methoas  give  results  which  are  in  very  close 
agreement  but  that  at  a  time  equal  to  two  tenths  of  a  second 
there  is  a  small  divergence.  This  small  divergence  is  due 
to  terms  of  oruer  higher  tiian  tne  fourth.  These  higher  order 
terms  must  be  taken  into  account  in  the  analytical  solution  for 
values  of  t  greater  than  about  0.10  second.  Wuen  t  is  smaller 
than  0.10  second  the  results  given  by  equations  (62),  (63), 
and  (04)  are  sufficiently  accurate.  Computation  ox  the 
coefficients  of  the  terms  higher  than  tne  fourth  order  would 
be  extremely  onerous.  The  best  procedure  is  to  carry  out  the 
numerical  integration  first  and  follow  with  a  check  by  the 
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analytical  solution  until  a  value  of  t  is  reached  such  that  terms 
of  order  higher  than  the  fourth  can  not  be  neglected.  If  the 
values  up  to  tiiat  point  as  obtained  by  the  two  methods  agree 
with  each  other  to  within  an  assigned  limit  of  divergence  the 
numerical  integration  can  be  considered  as  checked. 

For  the  commutation  of  amplitude,  the  aerodynamic  coefficients, 

Kjj  and  %  +  ^  ^  were  necessarily  assumed  as  no  data  were 

available  from  which  these  coefficients  could  be  calculated  for 
the  BOMB,  Demolition,  1000-lb.,  M44.  The  stability  test  v/ill 
subsequently  make  it  possible  to  calculate  accurate  values  of 
the  aerodynamic  coefficients  for  this  bomb.  However,  for  the 

present  study,  values  of  K-j  and  X.,  ♦  X^  ^  for  bomDs  of  similar 

shape  were  used  in  order  to  compute  the  approximate  cutting  point 
of  the  ring  necessary  to  produce  the  desired  amplitude  of  yaw. 

kg 

The  values  of  and  Kj-  +  for  the  BOMB,  Demolition  ,10CQ-lb. , 

M44  were  assumed  to  be  the  saute  as  those  for  the  M30.  The 
coefficients  for  the  latter  bomb  are  smaller  in  absolute  value 
than  the  corresponding  coeiiicients  for  other  bombs  of  the  current 
M  series. 


The  values  of  s  anu  n2  were  com out od  by: 

T 


Knd2  PQd* 


2s  =  kl  + 


Table  III  shows  the  arrangement  for  the  computation  of 
amplitude.  The  values  of  ana  eQ  are  obtained  from: 


€o  =  “  <Pt  " 


£0  =  “  - 


hcpT  cos  (pj  sin  q,T 


0 

h((^  cos  qv  -  qt2  sin  <pT)  -  l  C^V  sin  ^T+  cos  ^ 


A  is  computed  from: 


ft*  +  2usee  +  e2u2h2 
o  OOP  00 

u  n/  n2  -  s2 


(26) 


Table  III  was  continued  past  the  value  of  t  at  which  A 
became  ecuai  to  65°  or  1.1345  raciaiu.  Tne  value  of  t  for  wnioh 
A  d  came  equal  to  1.1345  rauians  was  calculated  by  inverse 
interpolation.  In  this  case,  =  0.233b  sec.  The  values 
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or  o>0  =  96.4861®  (66) 

The  value  o  ,  however,  does  riot  Include  the  effect  of  the 
o 

static  friction  between  the  nose  ^lug  and  the  retaining  ring  on 
the  suspension  cord.  Dr.  A.  G.  Charters  of  the  Ballistic  Research 
Laboratory  has  recommended  that  the  cut  end  oi  the  nose  plug  be 
rounaea  off.  The  paint  on  the  plug  was  reraovea  and  both  the 
suspension  ring  and  nose  plug  were  filea  and  then  smoothed  with 
No.  1  sane  paper.  Finally  both  tae  suspension  ring  ana  the  nose 
plug  were  thoroughly  lubricated  with  a  graphite  base  grease.  In 
this  manner,  the  static  friction  was  reduced  to  a  minimum.  Dr. 
Charters  also  suggested  bearing  the  nose  plug  on  a  roller  or 
ball  race  attaches  to  a  tain  roc.  instead  of  suspension  by  ring 
anu  cord.  This  wo  ale.  eliminate  the  static  friction  ana  only  the 


T 


effect  of  the  much  smaller  rolling  friction  would  have  to  be 
dealt  v/ith.  In  this  instance!  however,  it  is  estimated  that  the 
static  friction  will  cause  approximately  7°  delay  before  the 
ring  begins  to  slip  along  the  plug.  The  correct  value  for  the 
cut  angle,  o,  is  therefore  given  by 


CD 

corr. 
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JJhis'  report  describes  a  eurreirtr  method  for  obtaining  a  large 


initial  amplitude  of  yav/  for  tests  recently  conaucted^unuer  the 
direction  of  Colonel  L*  E.  Simon.  Tne  numerical  data  u or  the 


30)®!  Demolitionf  1000-lb.!  M^h^ve  been  furnished  tlie  writer 
by  ;nr.  E.  S.  Martin.  The  raatheuatxcui  portions  of  the  work 
have  been  verified  by  Mr.  Herbert  Roistacher  ana  Mr.  Simon 
Mowshowitz.  The  numerical  integration  oi  the  equations  of 
motion  lias  been  performed  by  Miss  Elsie  Magnus on. 


The  writer  is  indebted  to  Dr.  L.  S.  bederick,  Dr.  E.  J. 
McShane  and  Dr.  A.  C.  Charters  for  valuable  suggestions  and 
criticism. 


F.  V.  Reno 
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